Two-particle bosonic-fermionic quantum walk via integrated photonics 
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Quantum walk represents one of the most promising resources for the simulation of physical 
quantum systems, and has also emerged as an alternative to the standard circuit model for quantum 
computing. Here we investigate how the particle statistics, either bosonic or fermionic, influences a 
two-particle discrete quantum walk. Such experiment has been realized by exploiting polarization 
entanglement to simulate the bunching-antibunching feature of non interacting bosons and fermions. 
To this scope a novel three-dimensional geometry for the waveguide circuit is introduced, which 
allows accurate polarization independent behaviour, maintaining a remarkable control on both phase 
and balancement. 
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In the framework of quantum information processing, 
quantum walk has attracted much attention in the last 
few years [1]. Quantum walk is an extension of the clas- 
sical random walk: a walker on a lattice "jumping" be- 
tween different sites with a given probability. The fea- 
tures of the quantum walker are interference and super- 
position which lead to a non-classical dynamic evolution. 
Two different cases may be considered: discrete- and 
continuous-time quantum walks ^J. The properties of 
these two walks have shown several similarities how- 
ever, the discrete quantum walk exhibits a higher flex- 
ibility due to the possibility of tailoring the quantum 
coiriproperties to investigate different dynamic scenar- 
ios [3, 111 . By endowing the walker with quantum proper- 
ties, many new interesting effects appear: quantum walks 
allow the speed-up of search algorithms [6| and the real- 
ization of universal quantum computation Moreover 
it has been recently shown that quantum walks with a 
large number of sites exhibit a highly nontrivial dynam- 
ics, including localization and recurrence 0, HI- Within 
this scenario, a possible application is the investigation 
on biophysical systems, like the energy transfer process 
within photosynthesis [9]. 

Single-particle quantum walks yield an exponential 
computational gain with respect to classical random 
walks; it can be noted that they have an exact map- 
ping to classical wave phenomena and therefore they can 
be implemented using purely classical resources. On the 
other hand, quantum walks of more than one indistin- 
guishable particle can provide an additional computa- 
tional power that scales exponentially with the resources 
employed. This could be used to improve simulation per- 
formances in complex tasks, e.g. the graph isomorphism 
problem [10]. However, they need quantum resources to 
be implemented, since classical theory no longer provides 
a sufficient description. 

Different experimental implementations of single par- 
ticle quantum walks were performed with trapped atoms 
ions [i3, [Tsl, energy levels in NMR schemes Q, 



photons in waveguide structures [15], and in a fiber loop 
configuration 0, |5| . Very recently quantum walks of two 
identical photons have been performed [TgI, U^- How- 
ever, up to now no experimental demonstration on how 
the particle statistics, either bosonic or fermionic, influ- 
ences a two-particle quantum walk has been reported. 

In this work, we report on the implementation of a dis- 
crete quantum walk for entangled particles. By changing 
the symmetry of entanglement we can simulate the quan- 
tum dynamics of the walks of two particles with bosonic 
or fermionic statistic. These results are made possible 
by the adoption of novel geometries in integrated opti- 
cal circuits fabricated by femtosecond laser pulses, which 
preserve the indistinguishability of the two polarizations 
as well as provide high phase accuracy and stability. In 
the discrete quantum walk the walker is represented by a 
quantum particle -such as an electron, atom or photon- 
with an additional degree of freedom spanning a two- 
dimensional space and named the "quantum coin" (QC). 
At any given time the particle may be in a superposition 
of the two basis states, up (|[/)) or down {\D))^ represent- 
ing the two "coin faces" (Fig. [1^). The QC state directs 
the motion of the particle and the stochastic evolution 
by an unitary process. A key difference with the classi- 
cal case is that the many possible paths of the quantum 
walker may exhibit interference, leading to a very dif- 
ferent probability distribution of finding the walker at a 
given location. 

Indeed, the evolution of the walk can 
be described with the step operator E = 
Ej \j - 1)01 ® \U){U\ + \j + 1)01 ® \D){D\, where 
\j — and \j + l)(j| stand for the operators which 

move the particle in the higher and lower position of 
the lattice, respectively. The coherent action of the 
step operator E and coin tossing leads to entanglement 
between the position and the internal degree of freedom. 
After several steps, the counterintuitive profile of the 
quantum walk probability distribution emerges as a 
result of quantum interference among multiple paths. 
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FIG. 1: (a) Unidimensional quantum walk: depending on the 
result of the coin toss the walker moves upward (U) or down- 
ward (D). (b) Scheme of an array of beam splitters (BSs) for 
a four- steps quantum walk. Vertical dashed lines indicate the 
steps T of the quantum walk and the horizontal strips repre- 
sent the position \j) of the walker. In an array with an even 
(odd) number of steps the output ports J are grouped into 
the even (odd) final positions \j) of the walker, (c) Different 
behaviours of bosons and fermions on a BS. 



More complex distributions arise when two particles 
are injected into the same quantum walk. An important 
perspective is to combine in the same platform quantum 
walk and entangled states [HI, Ull . In this case both ini- 
tial and final state of the walkers are entangled states 
and, depending on the entanglement symmetry, differ- 
ent final distributions may be observed. By changing 
the symmetry of the entangled state we can simulate 
the quantum walk of two particles with integer or semi- 
integer spin [29]. The fermionic and bosonic behaviors 
drastically influence the dynamics of their quantum walk 
on the line. 

The simulation of single particle quantum walks on a 
line can be implemented using single photon states, beam 
splitters, phase shifters, and photodetectors [2l|,[22|. The 
quantum dynamics is achieved by an array of balanced 
beamsplitters (BSs) as shown in Fig. [TJd, each vertical 
line of beam splitters representing a step of the quan- 
tum walk. Horizontal strips represent the position \j) 
of the walker. If a photon, at time T and in the strip 
j, is incident downward \D) (upward \U)) on the BS we 
can represent its state as |j, 1^)t (b, ^)t)- The tran- 
sition from time T to time T + 1 is given by the BS 
operator \j,D)T ^ ^(|j - 1,1^)t+i - \j + 1,/7)t+i), 

1-^' ^ 72 d-^' + 1' + \^ - 1' ™s oper- 

ation simultaneously implements the coin (precisely the 
Hadamard coin C = (1 ^i)) and step operator E. 
Note that, if the particle starts at position \j = 0), at 
even (odd) times it will occupy only even (odd) positions. 

Provided that all the optical devices used in the walk 
are polarization insensitive, the polarization degree of 
freedom may be exploited to entangle the photons in- 
jected into the BS arrays. Moreover by changing the en- 
tangled state from a symmetric one, such as the triplet, 
into an antisymmetric one, the singlet, it is possible 
to mimic the quantum dynamics of two non-interacting 



bosonic and fermionic particles (see Fig. [Tb). It must be 
noticed that the experimental realization of such a net- 
work of BSs is exceedingly difficult with bulk optics, even 
for a small number of steps, since it requires a quadrat- 
ically growing number of elements. Furthermore, for a 
correct operation of the quantum walk, the phase intro- 
duced by the optical paths, passing from each beam split- 
ter to the following, must be controlled and stable. 
Our approach exploits an integrated waveguide architec- 
ture, which allows to concentrate a large number of opti- 
cal elements on a small chip and to achieve intrinsic phase 
stability due to the monolithic structure. In a waveguide 
implementation BSs are replaced by directional couplers 
{DCs) i.e. structures in which two waveguides, brought 
close together for a certain interaction length, couple by 
evanescent field. 

To realize the integrated optical circuits we adopted 
the femtosecond laser writing technology [H, . Briefly, 
nonhnear absorption of focused femtosecond pulses is ex- 
ploited to induce permanent and localized refractive in- 
dex increase in transparent materials. Waveguides are 
directly fabricated in the bulk of the substrate by trans- 
lation of the sample at constant velocity with respect 
to the laser beam, along the desired path. Since it is 
a single-step and maskless process, this technique allows 
rapid and cost-effective prototyping of new devices. Fur- 
thermore, it has intrinsic three-dimensional possibilities 
which have indeed been exploited in this work. 

In a previous work we have demonstrated that fem- 
tosecond laser technology can produce high-quality 
waveguides able to support polarization entangled pho- 
ton states [26|. In these femtosecond laser written 
waveguides, birefringence is low and does not affect sig- 
nificantly the coherence of the photons. Anyway the 
guided modes for the two polarizations are still slightly 
different (as shown in Fig. [2^) and this results in a resid- 
ual polarization dependence in the properties of the fab- 
ricated DCs. In fact, the coupling coefficient depends on 
the overlap integral between the two guided mode pro- 
files and is indeed quite sensitive to even small differences 
in the mode dimensions When several devices are 

cascaded, as in the case of a i^C array for implement- 
ing a quantum walk, small differences in the splitting 
ratios would accumulate and in the end affect the indis- 
tinguishability of the two polarizations. 

The three-dimensional capabilities of the present tech- 
nology can be exploited to tailor the polarization be- 
haviour of the DCs. We have fabricated several DCs with 
the waveguides lying on planes at different angles 6 with 
respect to the horizontal (see inset of Fig. [2)3), but fixed 
interaction length and spacing between the waveguides. 
It can be observed that the ratio between the measured 
coupling coefficients for the two polarizations {Cy and 
Ch) depends on 6>, as shown in Fig. [2)3. In particular 
there exists an angle for which the ratio between the two 
coefficients is one, i.e. the coupler becomes polarization 
insensitive. In order to realize an experimental imple- 
mentation of a discrete quantum walk with photons, we 
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FIG. 2: Integrated optical circuits, a) Measured inten- 
sity profile for the guided modes with polarization V and H, 
at 806nm wavelength. The 1/e^ dimensions are reported, h) 
Ratio of the estimated coupling coefficient for polarization V 
{Cv) and polarization H {Ch) in directional couplers fabri- 
cated with different angles between the waveguides (see in- 
set), but fixed interaction length (3mm) and distance (ll/xm). 
The fitting line is a guide to the eye. c) Schematic of the net- 
work of directional couplers fabricated for implementing a 4 
steps quantum walk. The color coding indicates the writing 
depth of the waveguides, which is varying from point to point. 
d) 3D representation of the basic cell of the network, which 
acts as a Mach-Zehnder interferometer 

fabricated a network of DCs^ all realized with the tilted 
geometry described above, where the two waveguides are 
brought at 11 /im distance, at an angle of 62°, in the in- 
teraction region, thus guaranteeing the polarization inde- 
pendence. The length of the interaction region is chosen 
as L = 2.1mm in order to obtain a balanced splitting 
ratio. 

In the interaction region the two waveguides are at 
different depths in the glass. To connect one coupler to 
the following, we designed a structure where the waveg- 
uides continuously vary the depth, as shown with the 
color codes in Fig. [2]3. The basic cell of the network, 
depicted in Fig. [2]i, acts as a Mach-Zehnder interferom- 
eter. For the correct operation of the quantum walk all 




FIG. 3: Experimental setup can be divided into three parts, 
(i) The source: polarization entangled photon pairs at wave- 
length A = 806nm generated via spontaneous parametric 
down conversion in a 1.5mm /3-barium borate crystal (BBO) 
cut for type-II non-collinear phase matching, pumped by a 
CW diode laser with power P — 50m [2^. Waveplates 
(WPs) allow generation of any single photon state and the 
Bell states A delay line is inserted to control the temporal 
superposition of the photons, which are injected into the in- 
tegrated device through single mode fibers (SMFs). Interfer- 
ence filters determine the photon bandwidth A A = 6nm. (ii) 
Integrated quantum walk circuit realized by ultrafast laser 
writing technique (see inset on the right), (iii) Measuring ap- 
paratus: the chip output is divided by a beam splitter (BS) 
and magnified through a set of two lenses. The photons cou- 
pled to multimode fibers (MMFs) are then delivered to single 
photon counting modules. The MMFs are mounted on mo- 
torized translation stages in order to select an arbitrary com- 
bination of two output ports and measure two-photon coinci- 
dences. Polarization controllers (PC) are used before the chip 
to compensate polarization rotations induced by the fibers. 



the interferometers present in the network must be phase 
balanced. This is intrinsically achieved with the highly 
symmetric three-dimensional geometry implemented in 
the network (Fig. [St and[2]i). The two central waveg- 
uides of the structure start with an initial separation of 
250/im to couple the device with a single-mode fiber ar- 
ray, while at the output the waveguides are separated by 
TOjam. The whole chip is 32mm long. 

To carry out and characterize the different quantum 
walks we adopted the experimental apparatus reported 
in Fig. [3l Different single photon and two photon 
states were injected in the network of DCs. The 
singlet-triplet transition within the Bell basis was per- 
formed by applying a 7r-shift in the phase (j) of the state 
^{\H)A\V)B^e''^\V)A\H)B) through rotations of H WP 
and QWP (see Fig. [3]). Fine phase adjustment was per- 
formed by a voltage liquid crystal (LC) device inserted 
on mode kA- 

The output of the integrated device is collected by a 
suitable telescope, splitted through a bulk beam splitter 
and then coupled to two multimode fibers (MMFs). By 
independently translating the MMFs on the arms C and 
D we select the output ports to detect, respectively / 
and J, and measure the single photon signals Sc{I) and 
Sd{J) and two-photon coincidences Ccd{I-, J)- 

As a first measurement we characterized the quan- 
tum walk circuit with single photons injected in either 
mode kA OT By measuring the output signals Sc{I) 
we obtained the single particle distributions. In or- 
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der to demonstrate the polarization insensitivity of our 
device, we repeated this measurement by injecting hght 
in different polarization states -horizontal, vertical, di- 
agonal and antidiagonal- always observing very simi- 
lar distributions [29]. Each experimental distribution 
was compared with the expected one by the similar- 

ity S = (E.,, ^Di,D[^^f/ Ai Eij which is 
a generalization of the classical fidelity between two dis- 
tributions D and D' . The mean value over the tested 
input polarization states is Siph = 0.992 ± 0.002. 

As a second step we injected two-photon entangled 
states. The distribution of the triplet and singlet states 
|^+) = 0) and |^~) {(f) = tt) emerging from the 
quantum walk was reconstructed by measuring the co- 
incidence counts Ccd{I-,J) for each combination of the 
indices / and J. The measured bosonic and fermionic 
distributions compared with expected ones are reported 
in Fig. [4^-b for |<l>+) and |^~). As done for the mea- 
surements on the sigle photon quantum walk, we plotted 
the probability distributions for the walkers to be in the 
final positions z, j of the quantum walk which is related 
to the probability of photons to emerge from the output 
ports / and J of the BS array p9j . 

In Fig. [Hwe observe how the particle symmetry affects 
the quantum walks. Note that some of the diagonal ele- 
ments of the fermionic two-particle walk are nonzero both 
in the theoretical and experimental distribution. Indeed 
Fig. m reports the probability distribution of the walk 
positions and not of the physical spatial modes. In fact, 
the expected probability to measure two fermions over 
the same output spatial mode is vanishing [29]. 

Furthermore, by considering a generic phase (j) (dif- 
ferent from and tt), it is possible to simulate the be- 
haviour of the quantum walk of two anyons, particles 
with a non semi-integer spin that represent a generaliza- 
tion of fermions and bosons [30]. Precisely, the entan- 
gled state |*^) = -^{\H)a\V)b + e''^\V) a\H) b) simu- 
lates two anyons characterized by creation operators sat- 
isfying CiCj = e^^CjCi and qc^ = e'^^CjCi + Sij. These 
systems exhibit both bunching and anti-bunching be- 
haviours (i.e. diagonal and off-diagonal elements in the 
final distribution). As a further measurement, we there- 
fore prepared some anyonic states 1^"^), in particular with 
= f , f , f TT, and measured the output probabilities. In 
Fig. St the distribution for = ^ is reported as an 
example of an anyonic behaviour. 

The experimental data can be compared with the the- 
oretical distributions by the similarity obtaining Sbos = 
0.982 ± 0.002 and Sfer = 0.973 ± 0.002 for the bosonic 



and fermionic quantum walk and 5, 



7r/4 



0.987 ±0.002, 



Sail = 0-988 ± 0.001 and Sani^ = 0.980 ± 0.002 for the 



anyonic quantum walks with 



f , f , fTT, respectively. 



The obtained results are in good agreement with the ex- 
pected behaviours. 

In conclusion, we presented the behavior of a discrete 
quantum walk based on an integrated array of symmet- 



ric, polarization insensitive, directional couplers in which 
Theory Experiment 
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FIG. 4: Two-particle quantum walks: ideal (left) and mea- 
sured (right) distributions of (a) bosonic, (b) fermionic and 
(c) anyonic (with = vr/2) quantum walks. 



two-photon polarization entangled states are injected. 
Exploiting the different statistics of singlet and triplet en- 
tangled states, such scheme allowed us to simulate how 
symmetric and antisymmetric particles travels through 
the quantum walk. 

The insensitivity to photon polarization, high-accuracy 
in the phase control and intrinsic scalability of the in- 
tegrated multi-DC network presented in this work, pave 
the way to further advanced investigations on complexity 
physics phenomena. For instance, by introducing suit- 
able static and dynamic disorder in the walk it would be 
possible to simulate the interruption of diffusion in a pe- 
riodic lattice, like Anderson localization I3ll-l34 . and the 
noise-assisted quantum transport effect |9l, [35j. 

During the final revision of this manuscript, the simula- 
tion of quantum statistics with entangled photons within 
a continuous quantum walk has been reported onhne ^ . 
This work was supported by EU-Project CHISTERA- 
QUASAR, FIRB-Futuro in Ricerca HYTEQ and PRIN 
2009. 
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Supplementary Information 

Chip fabrication and single photon characteri- 
zation - Waveguides were micromachined in borosilicate 
glass substrate (Corning EAGLE2000) using a cavity- 
dumped YbiKYW mode-locked oscillator, which delivers 
300/s pulses at IMHz repetition rate. For the waveg- 
uide fabrication pulses with 220nJ energy were focused 
170/im under the glass surface, using a 0.6N.A. micro- 
scope objective while the sample was translated at a 
constant speed of AOmm/s by high precision, three- axis 
air-bearing stages (Aerotech FiberGlide 3D). Measured 
propagation losses are 0.5dB/cm and coupling losses to 
single mode fiber arrays at the input facet are about 
IdB. Birefringence of these waveguides is on the order 
of B = 7 X 10~^, as characterized in f\\. 

In order to evaluate the effect of this residual birefrin- 
gence we measured the single photon distributions for 
horizontal, vertical, diagonal and antidiagonal polarized 
light. The obtained probability distributions for photons 
injected in both the modes kA and ks are reported in 
Fig. [5l It is clear from the plot that the behaviours 
are polarization independent, as confirmed by the simi- 
larities associated to each distribution, whose values are 
reported in table H both for photons injected in mode kA 
and kB- 



Input state 


Similarity kA 


Similarity ks 


H 


0.991 ±0.002 


0.992 zb 0.002 


V 


0.997 ± 0.002 


0.993 =b 0.002 


+ 


0.983 ± 0.002 


0.991 ±0.002 




0.995 ± 0.002 


0.992 ± 0.002 



TABLE I: Similarities between the experimental distributions 
and the theoretical one for horizontal (H), vertical (V), diag- 
onal (±) and antidiagonal (-) polarized input photons in the 
single-particle regime. 

Equivalence between entangled states and 
boson-fermion quantum walk . Let us consider a 
generic T step single-photon quantum walk implemented 
as in Fig. lb. If the particle is injected at mode 
I J), with input state aj|0), the walk performs a uni- 
tary transformation on the creation operator, namely 
ttj -> XIk^-^^^k- Let's now consider two photons in- 
jected into the walk in a polarization entangled state 
\\^jfj) = -^{alb'^j ± a'^jb\)\0) where at and 6t are the 



horizontal and vertical polarization creation operators, 
respectively The evolution of the walk becomes 



1^ 



ij/ 



1 



{UiKUjL±UjKUiL)a^j,b{\0) (1) 



The probability of detecting one photon at K and the 
other at L (without measuring the polarization) is given 
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FIG. 5: Measured output probability distributions and the- 
oretical expectations for a 4 steps quantum walk of a single 
particle injected into mode kA (a) and mode ks (b). The 
experimental results demonstrate that the DCs network leads 
to the same quantum dynamics independently from the initial 
polarization state and in agreement with the theory. 



by 



p±{I,J;K,L) = 



mj,KL 



fj,KL\^ 



for L^K 



for L = K 



(2) 



It is easy to show that J; L) and J; L) 

respectively correspond to the probabilities of detecting 
at positions K and L two bosons or fermions injected into 
the modes / ad J of the quantum walk. In fact, if at the 
input state we have two identical bosons characterized by 
commuting creation operator d\ and the output state 

is 44lO) ^ \ T.K,L^ij,KLd^KA\^)' The probability of 
detecting one boson at K and the other at L is precisely 
J; iiT, L). If at the input state we have two iden- 
tical fermions characterized by anticommuting creation 

llo) 



IJ,KL'-K'-L 



cLc 



operator c\ we have c|cj|0) \'}Zk 
and the probability of detecting one fermion at K and 
the other at L is p-{I^ J; L). Note that the dynam- 
ics of two-particle quantum walk cannot be reconstructed 
by multiplying the output probabilities {\Ujk\'^) of two 
single-particle walks. 

Furthermore it is evident from Fig. [5] that the plotted 
probability distributions refer to the position j of the 
walker and not to the output sites J of the BS array 
Indeed, referring to Fig. lb , for a walk with T* steps, the 
relation between the probabilities of photons emerging 
from one of the TV = 2T* outputs of the BS array {PP) 
and the final position of the walker is: 




FIG. 6: Theoretical (top) and measured (bottom) distributions for two-photon states emerging from the integrated array of 
BS. Depending on the symmetry of the two-photon Bell state injected into the device bosonic (a), fermionic (b) and anionic 
(c) behaviors are observed. 
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Then it is clear that in an array with an even (odd) num- 
ber of steps the output ports J are grouped into the even 
(odd) final positions j of the walker. 

These relations between the output probabilities of the 
BS array and the sites of the QW cause the diagonal 
elements of the two-particle fermionic distribution to be 
non- vanishing, as clearly shown in Fig. [G^b). This can 
be ascribed to the relation between the outputs Js of 
the BS array and the positions js of the walker: the 
symmetrization postulate of quantum mechanics makes 



the probability of the photons emerging from the same 
output of the BS array to be zero -i.e. Pff = 0-, all the 
same, according to the relations (j3j) generalized to the 
case of two-particle walk, a probability pj^^^^ ^ for the 
walk distribution can be observed. 

This is not surprising because two degrees of freedom 
are involved, the lattice site and the internal coin state: 
the antisymmetry for fermions in the same lattice site k is 
provided by the different internal coin state, one particle 
is in the \U) state and the other one in the \D) state with 
a total wavefunction of the form: 

|*fcfc) = -^[|j,C/)-U,D)] (4) 
which is clearly antisymmetric. 



[1] L. Sansoni et a/., Phys. Rev. Lett. 105, 200503 (2010). 



